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On basis of modification of Einstein's gravitational equations by adding the term f(R) cx (3R n , a geometric model 
of quintessence is proposed. The evolution equation for the scale factor a of the Universe is analyzed for the two 
parameters n = 2 and n = 4/3, which were preferred by previous studies of the early Universe. Another choice of 
parameters n and fj is proposed from the following reasons: the exponent n close to 1.2953 follows from the request 
for the evolution of the Universe after recombination to be close to the evolution of the flat FRW model with cold 
dark matter and reasonable age of the Universe defines the value of the coefficient (3 . Such a model corresponding 
to the evolution of the Universe with the dynamical A -term describes well enough the observational data. 



1. Introduction 

The discovery of accelerated expansion of the Universe 
[TJ [H [3] has stimulated the quest for mechanism of 
present inflation. The most famous theoretical model of 
dark energy (DE) is the cosmological constant A . The 
corresponding FRW solution for flat Universe with the 
present densities ratio for cold matter and dark energy 
(f2 m /fiA ~ 0.3/0.7) describes satisfactorily the evolu- 
tion of the Universe at low redshifts [21 HI [5]. However, 
the nature of the constant A -term has been remaining 
to be inexplicable during many years. 

It is well-known, the application of the constant A- 
term for modeling of the early Universe was initially 
confronted with principal difficulties. Solving the prob- 
lems of the very early Universe, this term has to be 
reducing by several order of magnitude during following 
evolution of the Universe. This problem was solved by 
rejection of A-term, and corresponding inflationary be- 
havior was determined by models of "effective A-term" - 
quasi-classical scalar fields in the one way. The progress 
of these models is well-known. Another way to describe 
the inflationary behavior is to take into account the po- 
larization of vacuum of quantum fields in the early Uni- 
verse. Taking into account of the effects of polarization 
leads to appearance of the terms non- linear on curvature 
in the Einstein-Hilbert action. In such models the infla- 
tion appears self-consistently. Let us note two issues: 

- a correction to the Einstein-Hilbert action with an 
arbitrary function of scalar curvature R is equivalent 
mathematically to the introduction of scalar field into 
the classical Friedmann cosmology [SJ [HI [H] ; 

- the terms of the form R n were investigated in the early 
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works done on the problem of singularity before obtain- 
ing of exact corrections to the Einstein action following 
from the one loop approximation [9 J . The part of such 
solutions approaches asymptotically to the Friedmann 
solutions with A-term, however physical results of the 
solutions have not been explained at that time. 

For the purpose to explain the accelerated expan- 
sion of the Universe today it is naturally to use the 
experience cumulative at investigations of the early Uni- 
verse. Thus, one of the tendencies is concerned with the 
hypothesis of existence some scalar fields that deter- 
mine the density of dark energy (see e.g. pH HU [12]). 
The another tendency models is an effective quasi- 
hydrodynamical energy-momentum tensor describing 
the observational data [13] , And the third tendency 
consists in generalization of the Einstein-Hilbert equa- 
tion by inclusion of curvature invariants [141 1151 I16j 
analogously to the earlier works mentioned above. The 
last approach one can consider to be either an inde- 
pendent approach to describing of DE or an analogue 
of inclusion of scalar fields (in case of f(R)) as stated 
above. 

In the works on the higher order gravity theories 
(HOGT), the models with power corrections were inves- 
tigated, however they have never been fitted to whole 
set of the observational data. 

In this paper, the model with correction f(R) oc R n 
with n > is considered in detail for the purpose to cor- 
relate it with the observational data. In the other words, 
we would like to obtain the model that does not conflict 
with the scenario of the large scale structure formation 
(in past) and describes satisfactorily the Universe under- 
going an accelerated expansion at present. Therefore, 
at the minimum, within the framework of the f(R)- 
theories, we will obtain solutions remind ACDM model 
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describing well enough by set of the observational data 
[3113]. However, as it was mentioned by various authors 
(see [5] and references therein), the observational data 
indicates the models of dynamical DE. Hence, our sec- 
ond aim is to search out such dynamical solutions within 
the framework of HOGT and to find out whether these 
solutions are preferable. 

This paper is organized as follows. In section 2 the 
basic equations of HOGT are presented. Section 3 is 
devoted to consideration of the models with n = 2 and 
n = 4/3. In section 4 for the corrections of the form 
f(R) oc (3R n we find the exponent n\ = 1.2953 which 
allows generalization of the Einstein equation for the 
scale factor a to have a particular solution correspond- 
ing to the flat FRW solution for cold matter. We show 
that instability of the solutions, that are close to this 
particular solution at z » 1, may lead to the acceler- 
ated behavior of the model at present and the following 
asymptotic approach of the solution to the solution with 
the constant A-term. In section 5 we discuss our results 
and compare them to observations. In the model there 
are only two free parameters - the coefficient (3 and a 
slight deviation of the parameter n from n\ = 1.2953 
mentioned above. Fixed from one set of the observa- 
tional data, they allow to obtain the rest of the set of 
the observational data. 

2. Basic equations 

As suggested by observations, we consider the flat cos- 
mological Friedmann model with the metric 



ds 2 = dr 2 - a{T) 2 (dx 2 + dy 2 + dz 2 ). 



(1) 



If Hq denotes the presently observable Hubble constant 
(the subscript will always indicate the present epoch), 
the reduced curvature tensor p\ = HQ 2 R L k has the fol- 
lowing matrix elements as a function of the reduced time 
9 = Hqt , with the notation a = da/d6: 



Pi 



— 3a/a, 

p = — 6 (a/a + a 2 /a 2 ) 



(2) 
(3) 



The variation of Einstein's Lagrangian with an addi- 
tional term AL(R) = f(R) gives 



(4) 



Here corresponds only to cold matter in the present 
Universe and 



rpk I ®f\ Jfc 1 rfc f 

= [aj)*-2 S *f- 



{5fg lm - 8\g km ) ( 



dp 



(5) 



specifies the effective quintessence with the nontrivial 
dependence on curvature. 



We proceed as in Ref. [T7J by introducing the new 
variable 



V = (aa) 2 



(6) 



which allows to reduce the order of the equations. Then 
the i = k = component of Eq. ([4]) leads to 



f, P v 



a dy 
2 da 



6 



f(p) + ay 



Pn 



dip 
da 



(7) 



fl m a, 



where p m is the a-dependent cold dark matter (CDM) 
energy density, p» = 3H 2 /8ttG is the critical density. 
Choosing the value of the scale factor a(#o) equal to 1 
at present, one has p m = po/a 3 and (p m //9*)a 4 = ri m a. 

In order to investigate the evolution of the cosmolo- 
gical model it is enough to obtain the solution of Eq. 
([7| with appropriate initial conditions. But for inter- 
pretation of the solution and for its comparison with 
observations it is necessary track for changes of CDM 
energy density p rn and quintessence energy density p v 
separately. For this purpose we use i = k = compo- 
nent of Einstein's equations j4| 



Lr 



1 



8ttG . 

-JpiPn 



Pv) 



(8) 



Multiplying this equation by a 4 and using Eqs. (HI),© 
and (0, we have 



V = {Pm + Pv)d 4 /p* 



(9) 



accounting for the evolution of cold matter energy den- 
sity. From this we find for the quintessence energy den- 
sity: 



— = (y - Q m a)/a 4 . 



(10) 



By solving Eq. ([7]) , we can find the evolution of p v . At 
known y, the Hubble parameter 

h(a) = y/jfi?, (11) 

the deceleration parameter of the Universe 
36 v w 



2 V 1 + 5. 
and the equation of state 
'2q- 1\ ( 5 V + 1 



5v = Pv/Pm 7 (12) 



(13) 

Here we will investigate corrections to the Einstein- 
Hilbert action of the form 

f(R) = -aR n . (14) 

In such a case, Eq. ([7]) can be presented in the form 



n{n-l)y"y + 



13 = (-3)' 



(1-n) 



(y'Y+n(4 
{y'?- n 



3n 



y y 



(y-Q m a), (15) 
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A general approach to investigation of the last equation 
is given in [T?] . 

From (|15p we see that a simple power in the asymp- 
totic solutions is absent for n — 2 and n — 4/3. These 
happen to be the same parameters which were of spe- 
cial importance in a previous theory of the early Uni- 
verse. At n = 4/3, one of us [5] has obtained for the 
first time a cosmological model without a singularity. 
That model passes through a regular minimum, has in- 
flationary stage and tends asymptotically to the classical 
Friedmann solutions. In addition, the usefulness of an 
additional term R 2 in the models of the early Universe 
has been pointed out before. Therefore we analyze the 
possibility of using of such powers for construction of 
models with variable parameters q and w . 



3. Models with 



n 



2 and 



n 



4/3 



These models have been considered in details in [To] , 
The model with n — 2 was often used in the theory of 
the early Universe. It describes the stage of fast oscil- 
lations of a(6) (the so called scalaron stage, which was 
introduced by A. Starobinsky in P~5]). The damping of 
such oscillations was connected with creation of unsta- 
ble particles and filling of the early Universe by a hot 
plasma. 

Here we want to consider this model in the opposite 
regime when the period of oscillations is commensurable 
with age of the Universe 1 / Hq . One can easily see the 
oscillations of the model by inserting the specified form 
of quintessence in (0} and ([5]) 



dO 2 



n adp 



0. 



(16) 



The scalar curvature performs oscillations near the value 
of p which corresponds to the model of cold dust matter 
in the Friedmann Universe. 

The special feature of the model with n = 4/3 con- 
sists in absence of the scale factor a in explicit form in 
the Einstein's equations if matter is neglected [9] . This 
allows one to find the general solution of Eq. (fT5j) with 
given n. One can show that de Sitter's solution arises in 
the limit a> 1. The curvature of this limiting solution 
is determined by the parameter a . 

For a further analysis of these models and compari- 
son with observable data see Section [37J 

4. The best-fit model 

In Eq. (|15j) there are two parameters, (3 and n, deter- 
mined by the observational data. These parameters 
could be chosen according to different requirements [151 
116] . Here we will choose parameter n from requirement 
of closeness of evolution of our model to the classical 
solution for the flat FRW Universe with cold matter in 
the past [T5]. This fact allows this scenario to be close 
to the scenario of the large scale structure formation. 



This requirement can be realized at condition that the 
classical Friedmann solution 



y 



a, 



(17) 



is a particular solution of Eq.(fT5j). It is easy to see from 
Eq.fT5j), the last condition is equivalent to the choice of 
n to be satisfying the equation 



ri-1 = 2n(4-3n) 
with roots 

7ii = 1-295, n 2 = -0.129. 



(18) 



(19) 



The first of the roots leads to the type of models of 
papers [15l [16] , while the second root corresponds the 
models with correction of the form oc /j,/R> n2 > investi- 
gated in [5D] . As it will be shown below, such a choice of 
n approaches the model to the set of the observational 
data in the best way. 

4.1. Behavior of the dust solution in the 
/(i?) -model in past 

After recombination, the evolution of the Universe has 
to be described by the Friedmann model with the cold 
matter. The dynamics of expansion is determined by 
stability of the dust solution in the model ([15]) . If the 
solution is stable, then the model evolves in the way 
very close to classical one. However, the observational 
data at z < 1 does not correspond to such a scenario, 
i.e. we are interested in dust-like solutions which are not 
stable in the model ([15]) but the perturbations do not 
grow catastrophically fast. Otherwise, the model does 
not provides a sufficiently long period with q ~ 0.5 at 
z 1 required for the large structure formation. 

For investigation of behavior of the solution ([TT]) we 
will search a perturbed solution in the form 



y = Qa(l + ip), ip<l. 



(20) 



In this linear approximation and close to the recombi- 
nation time (a <C 1) Eq. (fTS"]) yields, 

n a 2 ip" + aip'{2n- 0.5) = 0, (21) 

with the damping solution ip = C\ +C% a^" 1 -* . Hence 
the solution {TTJ) for the /(i?) -theory (TT3J) with n = Tii 
satisfying Eq. ([T8]) asymptotically approaches to to the 
flat FRW solution, i.e. it is stable and does not satisfies 
the requirement stated abov^l. 



3 It is interesting to note that the mentioned exponents (1181 arc 
obtained in the recent paper 121 1 for the equation equivalent to 
Eq.l ll5l l. In the /(i?) -theory, the given exponents allow obtain- 
ing of the solutions for cold dark matter coinciding with ones in 
the classical Friedmann model of the Universe. Also, it has been 
shown in 121 1 that these solutions are stable within the framework 
of HOGT. Let us notice we slightly change exponent n in the 
present work to obtain weakly unstable solutions adequate to the 
observational data. Authors are grateful to authors of paper 12 1 1 
kindly attracting our attention to their results. 
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As the next step on the path of the choice of n , we 
will look for solutions of the /(i?) -theories (fT4")l with n 
which is a little different from n\ , 

n = ni+Sn, 8n < 1. (22) 

It is efficient to rewrite the condition Q18p in the form 

S(n-1) = 2n(4-3n), 5 = 1 + e, e < 1. (23) 

We will look for the perturbed dust-like solution in 
the form (|2"0|) near to recombination. Eq. tfTir)]) for such a 
case yields 

n a 2 i>" + 2a ip' (n - 0.5 + 6/4) + {5 - 1) = 0. (24) 

Then after the change of variable £ = ln(a/a*) the Eq. 
(124]) yields 

nip + {n- 1/2 + e/2)ip + e = 0," = d/d£. (25) 
The last equation has a solution 

^( c '+t^K£D +c - ( *- i) - (26) 

The analysis of (|2"6"|) has shown that the requirement 
stated above is realized only for e < . 

In this case, the modification of the exponent (f2"3l) is 
determined by a small positive correction 

5n~ -2eni(m- l)/(12m- 1). (27) 

The numerical analysis has shown that behavior of so- 
lutions is sensitive to small changes of e at a — > 1 . The 
last fact together with a choice of parameter (3 allows 
us to obtain a good enough correspondence with the 
observational data. 

As an illustration we give the results for the set of 
parameters n = 1.296, (3 — 0.467 fixed according to 
Vt m = 0.29, h = 68 (the best fit to the CMB+SNe 
data presented in [S]). In a left panel of Fig. 1, a solid 
line represents the evolution of the variable y with the 
scale factor a. At the beginning, it coincides with the 
evolution of the dust model which is represented by a 
dotted line but further it deviates to the ACMD model 
represented by a dashed line. In a right panel, one can 
see evolution of the Hubble parameter with the scale 
factor a for the mentioned three models. 

4.2. The behavior of the solution in future 

The further expansion of the Universe at a > 1 ac- 
cording to Eq. (fT"5|) leads to the negligible effect of cold 
matter on the solution behavior. In this case, de Sitter 
solution y — il a 4 is an asymptotic solution of Eq. (fT5")) . 
This solution corresponds to the constant Hubble pa- 
rameter 

h(a -^oo) = VU (28) 
with £1 defined from equality 

3/3(2n- l)(n- l)/2 = (4fi)( 1 -"V4 (29) 



The inflationary solution is stable in the process of 
evolution of the model. To show it, we shall look for a 
solution with perturbation in the form 
y = fia 4 (l + <!>), $ < 1 . This ansatz yields 

n ia 2 $" -I- [(m - l)/2 + 3nf]o$' + (30) 
6(2n 2 - 3ni + 1)$ = 0. 

The change of variable a to variable £ (see Eq.(f26|) 
yields 

$ + A$ + B$ = 0, '=d/d£. (31) 

where coefficients are A = [3m — (1/ni + l)/2] = 2.99, 
B = 6(2m — 3 + 1/ni) = 2.18. This equation for per- 
turbations have a damped solution indicating de Sitter 
solution to be stable. The numerical analysis has shown 
that de Sitter solution is an attractive solution. 

As an example let us consider the evolution of the 
Hubble parameter for the case Q m = 0.29. In con- 
trast to the Hubble parameter of ACDM model mono- 
tonically decreasing down to constant v^Ta, it reaches 
a minimum h m « 0.978 at a w 1.15 and after that 
increases up to the asymptotic solution (|2"gj) . It is in- 
teresting to note that the formula for the dimensionless 
Hubble parameter h(z) = h{a) obtained from the obser- 
vational data in [5] allows its extrapolation to the future 
(a > 1) . At the parameters mentioned in this paper, the 
formula for h(a) also predicts minimum of the Hubble 
parameter at a w 1.45 which is equal to 0.951. The 
deceleration parameter q also passes a minimum and 
approaches to —1 with the growth of a. Therefore, we 
live in a transitional epoch between the classical Fried- 
mann cosmology and a de Sitter cosmology. 

5. Discussion 

Hereafter we will present the comparison of results of 
our model and the observational data. 

After fitting the model's parameters to the present 
observational data - acceleration of the Universe, the 
Hubble parameter at red-shift parameter z = Q, and 
the age of the Universe - there are no free parameters in 
the model. Its predictions for large z can be compared 
with observations. 

In case of n = 2 , the acceleration of the expansion 
changes at z ~ 0.46 to a deceleration (w w —0.72, q w 
0). Near z = 1 the variant of dust-like dark energy 
(w = 0) is realized. It corresponds to latest observati- 
onal data 22]. 

A special feature of the model is that for a — ► the 
variable y = (da) 2 tends to a constant which is equiva- 
lent to an evolution of the Universe filled by hot matter 
(p = e/3) in the Friedmann cosmology. When the so- 
lution is oscillating, the inflation at z ~ cannot be 
eternal although the period of oscillations is compara- 
ble with the age of the Universe. 
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In case of n = 4/3 , the crossover from a decelerated 
expansion to an accelerated one takes place at z ~ 0.30 
when w —0.60 . 

It is interesting to compare our model with the sim- 
plest ACDM model (for review see [22] ). Using the 
notation of section [23 we obtain 

V = On£ 1/2 + f« 2 (32) 

with the same values fi m ~ 0.3 and Q v ~ 0.7 

Our results show the following: the age of the Uni- 
verse is near 1/Hq « 13.7 Gyr. This age is larger than 
the age from the pure p 2 -model (« 11.4 Gyr) and from 
the p 4 / 3 -model (« 10.8 Gyr). The parameter w of the 
A CDM model calculating with use of ([3"3"| has the value 
— 1 . We may use Eq. (fT2"|) to determine the deceleration 
parameter q and compare it with the p 2 and p 4 / 3 mod- 
els. At z = 1, the deceleration parameter is q « 0.2, 
whereas the p 2 and p 4 / 3 models have q sa 0.5. This 
value is closer to the observational data. At z « 0.5, 
the parameter q in the p 2 -model is close to zero and in 
the ACDM model q sa —0.1. Apparently, the p 2 -modcl 
corresponds better to observations. The ACDM model 
is mathematically very simple. But it leaves the value of 
the A-term an unexplained fundamental constant. For 
dynamical models of the A-term (for example, the p 2 - 
model) this value evolves from a large Planck value in 
the early Universe to small value at present. 



As a set of observational data, the analysis of SNc 
and CMB data from [5] has been used. In that pa- 
per authors have reconstructed the resent history of 
the Universe on the base of SNe and CMB data in 
the model-independent way, only modeling DE by the 
hydrodynamical equation of state 

P = wp, w = (2q-l)/(3-n m /h 2 ). (33) 

The cited paper presents two conceptions of the analysis 
of the observational data: the first of them is the best fit 
to the data which uses only the hydrodynamical describ- 
ing of DE and does not impose restrictions on the values 
of il m and ho , while the second conception follows the 
priority of the concordance ACDM model, so authors 
of [5] put Q, m = 0.27 ± 0.04 and h = 0.71 ± 0.06. 

We will give the comparison of our results with both 
of them. Also we notice that the analogy of the "hydro- 
dynamical" DE ([33)) is not so proper to the higher order 
gravity theories, hence one can expect the comparison 
over the values h and q to be more informative than 
over the value w . 

As it has been found in [5], the best fit values are: 
Q m = 0.385, ho = 60. In the model DE evolves in time 
strongly enough. For given fl rn and ho we compared 
the results of the /(i?)-model with n = 1.2955 (Sn — 
0.0002) and (3 — 0.273 for the "geometric equation of 
state" parameter w = (y — ay')/3(?/ — fi m a) and the 
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Figure 2: The comparison of results of f(R) oc R n model with n = 1.296 and (3 — 0.467 (the thick solid line) to results of 
analysis of SNe+CMB data with Q m — 0.29 done in [5]. The evolution of the deceleration parameter with redshift is shown 
in a right panel and variation of equation of state of DE is shown in a left panel. The best fit of SNe+CMB data in such case 
is represented by the thin solid line, the la and 2a confidence levels are represented by the light and dark grey contours, 
respectively, and ACDM is represented by the dashed line. 



deceleration parameter q with the results of [5] . In this 
/(i?) -model the age of the Universe is 14.9 Hyr, the 
deceleration parameter is qo = —0.91 at present and 
the transition to acceleration occurs at z = 0.38. Simi- 
larly to results of [5], the wr>E < — 1 at lower redshifts 
(wdeo — —1-53), however , the evolution of equation of 
state of "geometric DE" is more weak contrary to the 
results of [5]. 

However, if strong priors have been imposed on n m 
and h (i.e. the ACDM model priors: fl m = 0.27±0.04 
and ho — 0.71 ± 0.06), the evolution of DE is extremely 
weak and in good agreement with the ACDM model. 
The best fit in the case is fi m = 0.29 and there is a 
good enough coincidence of our model and their analysis 
for parameters of the model n = 1.296 (Sn = 0.001) and 
(3 = 0.467 (see Fig. 2). The deceleration parameter qo = 
—0.683 at present, and the deceleration was changed by 
the acceleration at z — 0.51 (qo = —0.63 ± 0.12 and 
z = 0.57 ± 0.07 in |5J). The age of the Universe in this 
case is 13.6 Hyr. 

Thus, the f(R) oc j3R n -model with parameters (3 
and n chosen according to the principles mentioned 
in Introduction describes the evolution of the Universe 
quite corresponding to the SNe+CMB data. 
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